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ABSTRACT 

Binaries are not always neatly aligned. Previous observations of the DI Her system showed that the 
spin axes of both stars are highly inclined with respect to one another and the orbital axis. Here we 
report on a measurement of the spin-axis orientation of the primary star of the NY Cep system, which 
is similar to DIHer in many respects: it features two young early-type stars (~ 6Myr, B0.5V+B2V), in 
an eccentric and relatively long-period orbit (e = 0.48, P = 15^3). The sky projections of the rotation 
vector and the spin vector are well-aligned (/3p = 2 ± 4°), in strong contrast to DIHer. Although 
no convincing explanation has yet been given for the misalignment of DIHer, our results show that 
the phenomenon is not universal, and that a successful theory will need to account for the different 
outcome in the case of NY Cep. 

Subject headings: techniques: spectroscopic — stars: kinematics and dynamics — stars: early-type — 
stars: rotation — stars: formation — binaries: eclipsing — stars: individual (NY Cep) 



1. INTRODUCTION 

We are conducting measurements of the relative orien- 
tations of the rotational and orbital axes in close binary 
star systems. Our name for this undertaking is the Ba- 
nana project, an acronym chosen to remind us that bi- 
naries are not always neatly aligned. This paper reports 
our results for the NYCephei system, the third system 
we have studied , the f irst two having be en V1143Cy^ni 
(lAlbrecht et al.||2QQ7[) and DIHerculis ( |Albrecht et al. 
2UU9[ ). Those two earlier works were written before our 



project was enlarged and named, and are retrospectively 
included in this series as Papers I and II. The goal of 
the Banana project is to enlarge the number of detached 
binary systems for which the relative orientation of the 
spin axes is known, in order to shed light on the formation 
and evolution of binaries and perhaps also of planets. 

1.1. Motivation 

Close binaries and star-planet systems might be ex- 
pected to have well-aligned orbital and spin angular mo- 
menta, since all of the components trace back to the same 
portion of a molecular cloud. However, good alignment 
is not guaranteed. If a cloud is highly elongated, with its 
long axis tilted with respect to its rotation axis, then the 
cloud may give birt h to binary stars wit h a strong spin- 
orbit misalignment ( Bonnell et al.|p!992 ). Alternatively, 
disks around young stars might become warped during 
the last stage of accretion. This warp could torque the 
orbit by a large angle while m aintaining the orientation 



of the spins ( [Tremaine 1991). More generally, star for 
mation may be a chaotic process, wi th accretion from 



different directions at different times ( Bate et al.||2QlQ ), 



pbreclit@space.niit.edu 

Based on observations made with Sophie, a high-resolution 
echelle spectrograph on the 1.93-m telescope of the Observatoire 
de Haute-Provence. 



and perhaps we should not expect the angular momenta 
of the star and the disk to be as well aligned as they 
apparently were in the Solar system. 

There are also processes that could alter the stellar 
and orbital spin directions after their formation. A 
third body orbiting a close pair on a highly inclined or- 
bit can introduce large oscillations in the orbital incli - 
nation and eccentricity of the close pair (Kozai 1962). 
Tidal dissipation during the high-eccentricity phases can 
cause the system to free itself of these "Kozai oscilla- 
tions" and become stuc k in a high-obliquity state ( |Fab^ 
rycky fc Tremaine||2QQ7 ). However, if dissipation is suffi- 



ciently strong then the system will evolve into the double- 
sync hronous sta te, characterized by spin-orbit alignment 
(e.g. |Hut||1981[). Therefore, whether a close binary or 
a star-planet system is well-aligned or misaligned de- 
pends on its particular history of formation and evolu- 
tion. Even though this issue is important for a complete 
understanding of star formation, there has been very lit- 
tle observational input. 

Another motivation comes from exoplanetary science. 
Recently it was revealed that many of the close-in giant 
planets ("hot Jupiters") have orbits that are strongly 
misaligne d with the ro tation a xes of their parent stars 
see, e g.,[H6brard et al . 2008; Winn et al.||2009[ [Narita 



et al.| |2009[ |'l'riauT^al.| |201Qt . The high obliquities 



are especially common am ong stars with higher masses 
and effective temperatures (Winn et al. 2010 'Schlaufman] 
2010). The results of such studies are commonly inter- 
preted as constraints on theories of the "migration" pro- 
cesses that presumably brought the planets inward from 
their more distant birthplaces. However, some theories 
invoke processes that produce large stellar obliquities for 
reasons having no thing to do with t he planets, such as 
chaotic accretion (Bate et al. 2010), and magnetic in- 



teractions with the inner edge ot the accretion disk (Lai 
et al.|201Q ). Although those theories have not been fully 
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TABLE 1 

Measurements of spin orbit angles in eclipsing binaries 



System 


Spectral Type 


Period [d] 




Eccentricity (e) 






reference 


^Lyr 


Be+B6-8II 


12.9 


0.52 


<0.01 




aligned? 


1,2 


VlOlOOph 


A7IV-V+? 


0.66 


0.4511±0.011 


= 




aligned? 


3,4 


WUmi 


A3V+G9IV 


1.7 


0.363±0.001 


= 




aligned? 


5 


(5 Lib 


AOV+KOIV 


2.33 


0.30±0.06 


= 




aligned? 


6,7 


V505 Sgr 


A2V+GIV 


1.2 


0.288±0.001 


= 




aligned? 


8,9 


AIDra 


A0V+F9.5V 


1.2 


0.284±0.004 


= 




aligned? 


8,10 


XTri 


A3V+G3IV 


0.97 


0.272±0.003 


= 




aligned? 


5 


RZ Cas 


A3V+KIV 


1.20 


0.233±0.001 


= 




aligned? 


5 


USge 


B8.5V+G3III 


3.38 


0.219±0.002 


0.04 




aligned? 


5 


WW Cyg 


B8V+G4III 


3.32 


0.215±0.002 


= 




aligned? 


5 


YLeo 


A3V+K3IV 


1.69 


0.213±0.001 


= 




aligned? 


5 


RXHya 


A8+K0IV 


2.28 


0.211±0.013 


= 




aligned? 


5 


Algol 


B8+K2IV 


2.87 


0.206±0.003 


= 




aligned? 


11,12,13 


RWGem 


B5-B6V+F0III 


2.87 


0.198±0.002 


= 




aligned? 


5 


YPsc 


A3V+K2IV 


3.77 


0.193±0.001 


0.12 




aligned? 


5 


TV Cas 


A2V+G1IV 


1.81 


0.188±0.15 


= 




aligned? 


5 


ST Per 


A3V+KIV 


2.65 


0.184±0.004 


= 




aligned? 


5 


UCep 


B7V+G8III 


2.49 


0.177±0.009 


= 




aligned? 


5 


TXUma 


B8V+F7-F8III 


3.06 


0.164±0.002 


= 




aligned? 


5 


WDel 


A0-B9.5V+K0IV 


4.81 


0.151±0.004 


0.20 




aligned? 


5 




saw-D-t- ; 


U.ZD 


n 1 /L-i-n ni 
u. ±4tzizu.u± 


= u 




aligned? 


1 A 


DEDra 


B0V+B2V 


5.3 


0.14 


0.018±0.011 




misaligned? 


15 


SWSyg 


A2V+K1IV 


4.57 


0.138±0.002 


0.30 




aligned? 


5 


RYPer 


B4V+F0III 


6.86 


0.137±0.003 


0.21 




aligned? 


5 


RZ Set 


B2II+A0n-III 


15.19 


0.136±0.007 


= 




aligned? 


5 


AQ Peg 


A2V+K1IV 


5.55 


0.123±0.002 


0.24 




aligned? 


5 


RYGem 


A2V+K0III-K1IV 


9.30 


0.097±0.014 


0.16 




aligned? 


5 


NYCep 


B0V+B2V 


15.3 


0.086±0.015 


0.445±0.004 




/3p = 2° ± 4° 


16,17 


DIHer 


B4V+B5V 


10.55 


0.0621±0.001 


0.489±0.003 


^P 


= 72° ±4°, /3s = -84° ±8° 


18,19 


V1143 Cyg 


F5V+F5V 


7.64 


0.059±0.001 


0.5378±0.0003 




= 0?3±1?5, /3s = -1?2±1?6 


20,21 



Note. — In column 4, rp denotes the radius of the primary star divided by the length of the orbital semimajor axis. For DE Dra and 
(3 Lyr no uncertainties in the primary radii are given in the references. In column 5, the entry "= 0" indicates that the eccentricity was 
assumed to be zero by the authors. In column 6, /5pand /3s denote angles between the projections of stellar and orbital spin axes for the 
primary and secondary, respectively, using the coordinate system of Hosokawa (1953). 

References. — (l)^ossiter (1924); (2) Harmanec (2002); (3) Worek et al. (1988) (4) Corcoran et al. (1991); (5) Twigg (1979); (6) 



Bakis et al. 



( | 200 6|>; (^7 )]worek (1985); (8) Worek (1996); (9) Lazaro et al. (2UU6) 



vey (1931); (13j ISoderhielm (1980); (14) Kucinski (2009); (15) Mube . 



study; (IHj ^Popper (1982) ; (l9j Albrecht et al . (2009); (20j Andersen et a l. (1987); (21j Albrecht et al. (2[){)7j 



10) Lazaro et al. (2004); ( 1 1) McLaughlm (1924); (12) 
Couch (1982); (16) Mol mgrerTet al | ( |199D| rr (17) This 



developed, it would seem that the mechanisms they pro- 
pose should also operate in the case of binary stars, and 
therefore the theories might be tested by measuring the 
obliquities of binary stars. 

For example if interactions with a distant companion 
are important for close double stars, then this should also 
be true for star-planet systems. Other effects like tidal 
realignment do depend on the mass and mass ratio of the 
close pair. For a direct comparison with the exoplanet 
hosts, one would want to survey main sequence F and G 
binaries, while (as we will describe) most of the existing 
data, including the data presented in this paper, is for 
earlier type stars. We hope to rectify this situation with 
future observations. 

1.2. Measuring stellar obliquities 

Measuring the orientations of stellar spin axes is not 
straightforward. Telescopes cannot usually resolve stel- 
lar disks, causing the information on spatial orientation 
to be lost. An optical interferometer, in combination 
with a high resolution spectrograph, might allow for spa- 
tial resolution of rotationally-broadened stellar absorp- 



tion lines f Petr ov|1989l IChelh fc Petrov|1995] ) , but so far 



and brightest systems (Le Bouquin et al. 



such measurements are possible for only the ver y nearest 



2009). 



One would be able to determine the inclinations of the 
spin axes with respect to the sky plane, based on em- 
pirical estimates of the projected stellar rotation speed 
{vsini)^ the stellar radius (i^^), and the stellar rotation 
period (Prot), using the equation 



I = sm 



V smi 



(2^i?,/Prot) 



(1) 



Many investigators have pursued this path a nd found it 
to be blocked, for reasons described clearly by |Soderblom 
( |1985| . Namely, the quantities vsini, R^, and Prot 
are difficult to measure with high enough accuracy, and 
the flattening of the sine function near 90° prevents 
the method from discriminating even modest inclina- 
tion s from edge-on inclinations. D espite these difficul- 



ties, 



Glebocki & Stawikowski (1997) used this method to 



argue that active binaries that are observed to be asyn- 
chronously rotating have misaligned spin axes. 

Given a large sample of vs'mi measurements of stars 
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TABLE 2 
General data on NY Cephei 



HIP 
HD 

R-A.J2000 

Dec.j2ooo 

Vmax 

Sp. Type 

Period 

Eccentricity 

Inclination 

Rp 

Rs 

Mp 

Ms 

Ls/Lp 



tData from |ESA|fl997| > 
fData from Ahn (19 92|> 



113461 

217312 

22^58^40^ 

63°04'38'' 

7.5 mag 

B0.5V+ B2V 

15^27 

0.48(2) 

78(1)° 

6.8(7)i?0 

5.4(5)i?0 

13(1)M0 

9(1)M0 

0.37(6) 



:*rData from Holmgren et al 
Note. 



( 1990t 



Rp denotes the radius of the primary component 
and Rs the radius of the secondary component. Lp/Ls denotes the 
luminosity ratio between the primary and secondary. 




Fig. 1. — Pole-on view of the orbit of NY Cep. The solid 
and dashed ellipses indicate the orbits of the primary and secondary 
stars, respectively. The filled (red) circles indicate the positions of 
primary star and the open (blue) circles indicate the position of the 
secondary star during times of observations. The lines connecting 
the orbits with the center of gravity indicate the position of peri- 
astron. We observe the system from the x origin and a negative y 
value. 



that have a known distribution of rotation speeds, or 
that are presumed to have similar rotation speeds (e.g., 
main-sequence stars of a given mass and age), it is pos- 
sible to test for departures from an i sotropic d istribution 
of spin dire ctions. Using this idea, Guthrie| ([l985) and 



Abt (2001) searched for, and did not tind, a tendency 
for stars to be preferentially aligned with the Galactic 
plane. When the orbital orientation is also known, as is 
the case for visual binaries or eclipsing binaries, then such 
tests give information about spin-orbit alignment. Such 
studies are worthwhile but they are hampered by the 
nonlinearity of the sine function (as mentioned above) as 
well as uncertainty in the underlying distribution of rota- 
tional speeds, and the heterogeneity in the techniques for 
measurin g sin i. Thi s meth od has been used in various 



guises by|Weis|(|1974D,[Hale|(ri9 94J, and|Howe fc C larke 



( 2009 ) , among others, kecentlyj" JSchlaufman ( |2010 ) used 



this niethod on stars with transiting planets, tinding ev- 
idence that more massive stars have high obliquities. 

Indirect evidence on the orientation of young stars 
comes from the orientation of their disks - assuming star- 
disk alignment. The projected rotation angle of the disk 
can be traced by the linear polarization vect or of the star 
light reflected from the dust grains in disks (iMonin et ah 



1998). These aut hors and o thers (e.g. "Donar et al.'T999' 
Wolf et al. 200T| [Jensen e t al. 2004; Monin et al. 2006) 
mid that in most binary systems the disks around the in- 
dividual stars are aligned with each other. However they 
also mention exceptions to this rule. One of these excep- 
tions is the hierarch ical triple system T Taur i, for which 



Skemer etaTI (| 2008[ ) and [Ratzkaet al.| ( |2009[ ) found that 
c around the northern component, the original 



the disk 



T Tauri system, is viewed face on, while the disk around 
one of the southern components (TTauSa), separated 
by > 100 AU from T Tau N, is viewed edge-on. Thus this 
system has misaligned circumstellar disks. 

The method that provides the most accurate informa- 
tion for individual systems takes advantage of eclipses. 



During an eclipse of one star by another star or a planet, 
part of the rotating stellar surface is hidden, causing a 
weakening of the corresponding velocity component of 
the stellar absorption lines. Modeling of this spectral 
distortion reveals the relative orientation of the spin and 
orbital axes on the sky: the projected obliquity. This "ro- 
tation anomaly" was first predicted by Holtj (|1893|). A 



cl aim of its detection in the J Librae system was made 
by Schlesinger ( 1910|), but more definit ive measurements 



were achieved b'y'^Rossiter 
for the /3Lyrae and Algo 



1924j) and M cLaughlin| ( |1924 ) 
. systems, respectively. I'he 



phenomenon is now known as the Rossiter-McLaughlin 
(RM) effect. Various aspects of the theory of the ef - 



fect have been worked out by [Struve fc Elvey (|1931[); 
[Hosokawa (T953^; "Kopal' (T9590 ('1974[); |0'ht_a em. 

TO05); Gimenez (2006); Hadrava (200^ and |h 
[eFl ' ^ 



The work described in this paper. 



irano 
as well 

as in the previous papers in this series, is based on the 
RM effect. 

1.3. Previous observations of the RM effect 

Although it has been more than 80 years since its dis- 
covery, there are relatively few quantitative analyses of 
the RM effect observed in stellar binaries. It is ironic that 
at present, there are more such papers about exoplane- 
tary systems than about stellar binaries. In the past, 
observing the RM effect was generally either avoided (as 
a hindrance to measuring accurate spectroscopic orbits) 
or used to estimate stellar rotation speeds. Almost all 
authors explicitly or implicitly assumed that the orbital 
and stellar spins were aligned. 

Table [T] shows the results our literature search for anal- 
yses of the RM effect in stellar binaries. (This table also 
gives the results for NY Cep, the subject of this paper.) 
We have tried to be comprehensive, but cannot claim our 
list to be complete. The authors would appreciate being 
notified of any omissions. The systems are listed in order 
of decreasing rp, the radius of the primary star in units 
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OOOOX Z<2 OOOOX Z<2 



4440 4450 4460 4470 4480 4490 4500 
Wavelength [A] 

Fig. 2. — Example spectrum of NY Cep. Displayed is order 
number 11 of one spectrum of NY Cep, as delivered by the instru- 
ment software. No flat flelding, binning, or any other modiflcations 
to the data have been applied. The strong feature at 4471 A is the 
He I line. Adjacent to this line, at 4481 A, is the Mgll line. 



of the orbital semimajor axis, and thus the ordering is 
approximately a progression from closely-interacting sys- 
tems to well-detached systems (although this is not true 
for some of the Algol systems, which have large faint sec- 
ondaries). Column 6 summarizes the results for the stel- 
lar obliquities. In most cases we have written "aligned?" 
because the RM data appear visually to display the pat- 
tern of a well-aligned system — a redshift during the first 
half of the eclipse, followed by a blueshift of equal am- 
plitude during the second half of the eclipse — with the 
question mark indicating that no quantitative analysis 
was undertaken, and hence the uncertainty is unknown. 
For systems with orbital inclinations very close to 90° 
the results are especially ambiguous because in such cases 
there is a strong degenerac y between the projected obliq- 
uities and rotation rates ( Gaudi fc Winn |2QQ7[) . For 
DEDra we wrote "misalignedY" because |Hube ^ Couch 
(|1982") found the RM effect to be asymmetric about the 
mideclipse time, but gave no quantitative result for the 
projected obliquity. The only cases where quantitative 
results for the projected obliquities are given are from 
the Banana project. 



1.4. The Banana project and NY Cep 

For our program, we have begun by concentrating on 
relatively young, detached systems with < 0.15, in 
order to minimize the effects of tidal interaction, and 
thereby study a more "primordial" distribution of obliq- 
uities. We also omit Algol-type systems, to avoid the 
extra complexity of the spin evolution caused by mass 
transfer between the stars. These criteria exclude the 



Algol-type systems studied by Twigg (1979) as well as 
all the systems listed in Table [T] except tor DEDra, 
V1143 Cyg, DIHer, and NY Cep. The result of misahgn- 
ment in the DEDra system by Hube & Couch (1982) 
is intriguing but needs to be checked, l^br V1143 Cygni 
we showed that both s pin axes are wen-aligne d with the 
orbital axis [Paper I; Albrecht et al. (2007)]. In con- 
trast, for our second target, DIHer, we measured spin 
axes that are drastically misaligned with the orbital axis 
[Paper II; , Albrecht et al.,(2009j ]. This was the first clear 




4465 4470 4475^ 4480 
wavelength [A] 



4465 4470 4475^ 4480 
wavelength [A] 



Fig. 3. — Disentangled spectra of the two stars in the 
NY Cep system. The two upper panels show spectra between 
4461 A and 4484 A. The upper left panel is based on the sum of 
all the spectra obtained outside of the eclipse, after shifting into 
the rest frame of the primary star and after subtracting our best- 
fitting model of the secondary spectrum. The data are shown as 
gray lines, and the model as a continuous curve. The wavelengths 
of the lines used to construct the model spectrum are indicated 
above the panel. The lower left panel shows the difference between 
the data and the model. The right panels show the secondary 
spectrum, after subtracting the best-fitting model of the primary 
spectrum. For examples of the observed spectra with both stars 
present, see figures |6] and |7| 



demonstration of such a strong misalignment in a close 
binary. It resolved the longstanding problem of the sys- 
tem's anomalous apsidal motion. 

In this paper we focus on the NY Cephei system, whose 
properties are summarized in Table [2] A summary of the 
history of observations of this system was given by "AhnJ 
(1992). It harbors two early B type stars on an eccentric 
15^27 orbit. One peculiar characteristic of this system is 
the lack of secondary eclipses. Inferior conjunction oc- 
curs when the stars are relatively far apart, and with an 
orbital inclination of <80° the sky-projected separation 
of the stars is too large for eclipses. The orientation of 
the orbit is illustrated in Figure [l] 

Section [2] of this paper presents our observations and 
describes our model for the spectroscopic data, and our 
analysis procedure. The results for the orbital and stel- 
lar parameters, including the measurement of the stellar 
orientation are presented in section [3] The results are 
discussed in section lU 

2. OBSERVATIONS AND DATA REDUCTION 

We observed NY Cep with Sophie^ a high-resolution 
echelle spectrograph on the 1.93 -m telescope of the Ob - 
servatoire de Haute-Provence (iPerruchot et al.||2008|), 



employing its High Efficiency (k ~ 4U,UUUj mode! We 
chose an integration time of 20min resulting in signal- 
to-noise ratios (SNR) between 50 and 150 per pixel at 
wavelengths near 4, 500 A, the spectral region used for 
our analysis. To illustrate the data quality. Figure [2] dis- 
plays the region of interest of one of the spectra. 

The ffist half of the primary eclipse was observed dur- 
ing the night of 2009 Sep. 11/12. On that night, clouds 
prevented us from obtaining pre-ingress observations. 
The second half of the primary eclipse, as well as a num- 
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Fig. 4. — Illustration of the effects of rotation, macro- 
turbulence, pressure broadening, and eclipse blockage on 
the absorption line model. In each panel, the sohd hne rep- 
resents an absorption hne out of echpse, while the (red) dashed 
line shows the same line just before mid-eclipse. The left panel 
shows only the effect of stellar rotation. The middle panel shows 
broadening by rotation and macro-turbulence (as appropriate for 
the Mgll line). The right panel shows a line additionally affected 
by pressure broadening (as appropriate for the He I line). 



ber of spectra directly after egress, were observed during 
the night of 2009 Oct. 12/13. Strong wind and poor see- 
ing caused lower SNRs than we had anticipated. In addi- 
tion, we gathered some spectra on days spread through- 
out September to November 2009, in order to establish 
the spectroscopic orbit. All together, 27 spectra were 
obtained during primary eclipse, 8 spectra were obtained 
directly after the end of egress, and another 11 observa- 
tions at various orbital phases outside of eclipses, for a 
total of 46 spectra. We also attempted photometry of 
the eclipse throughout autumn and winter of 2009/2010, 
but all our attempts were foiled by bad weather. 

The two-dimensional reduced spectra were examinee^ 
and bad pixels were flagged and omitted from subsequent 
analysis. The spectra were shifted in wavelength to ac- 
count for the radial-velocity of the observatory relative 
to NYCep. Initial flat fielding was performed using the 
nightly blaze functions. We adopted the wavelength so- 
lution delivered by the spectrograph software, which has 
an uncertainty of order 1 ms~^ and is negligible for our 
purposes. 

2.1. Description of the model 

Our analysis focused on the 11th echelle order of the 
Sophie CCD, which encompasses the wavelength range 
from 4459 A to 4486 A, including the two best absorp- 
tion lines available for this study. He I (4471 A) and Mg II 
(4481 A). The helium line is strong, with a line width 
dominated by pressure broadening. The magnesium line 
is relatively weak but has the virtue of being chiefly 
broadened by rotation and therefore well-suited to the 
analysis of the RM effect. There are also a number of 
other weaker lines in this order that must be modeled si- 
multaneously with the stronger lines. The lines are illus- 
trated in Figure |3) which shows the spectrum of each star 
individually, after disentangling them with the modeling 
procedure described in this section. The spectrum from 
that order was binned to a resolution of about 9 kms~^, 
giving 213 pixels in the region of interest. The SNR 
was estimated from the scatter in the continuum on both 
sides of the He I (4471 A) and Mgll (4481 A) lines. 

^ The 2D-spectra can be obtained from the following web ad- 
dress: http : / / atlas . obs-hp . f r/ sophie/ 
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Fig. 5. — Results for the timing parameters, based o n our 
Monte Carlo analysis of the times of minimum light given by | Ahn| 
( |1992| ). The gray shading indicates the density of Monte Carlo 
results. The marginalized probability distributions are shown on 
the sides. 



Our spectral model is similar to the models described 
in Papers I and II. Because of the light from the fore- 
ground star, we cannot treat the RM effect as a simple 
wavelength shift, as is commonly done for eclipses of stars 
by planets. Instead we modeled the line profiles of both 
stars simultaneously, taking into account stellar rotation, 
surface velocity fields, orbital motion, and partial block- 
age during eclipses. By adjusting the parameters of the 
model to fit the observed spectra, we derived estimates 
of the orbital and stellar parameters. 

For each star, and for each phase of the eclipse, we 
created a discretized stellar disk with about 200,000 pix- 
els in a cartesian coordinate system. We assumed the 
disk to be circular, since the stars are well detached and 
have slow rotation speeds relative to the breakup veloc- 
ity. Each pixel has its own emergent spectrum, weighted 
in intensity according to a linear limb-darkening law, and 
Doppler shifted due the combined effects of orbital mo- 
tion, rotation, and macroturbulence. The orbital motion 
is specified by a Keplerian model common to all pixels. 
The rotation is assu med to be un iform (no differential 
rotation). Following Gray (2005), the macroturbulent 
velocity field is assumed to obey Gaussian distributions 
for the tangential and radial components, with equal am- 
plitudes, brightnesses, and surface fractions. 

The sky-plane coordinates of the foreground and back- 
ground stars are calculated, and if the background star 
is being eclipsed then the eclipsed pixels are assigned 
zero intensity. The emergent spectra from the uncov- 
ered pixels of both stars are summed to create a model 
absorption line kernel. For the He I line a further convo- 
lution with a Lorentzian function is applied, to account 
for pressure broadening. Figure |4] illustrates the effects 
of rotation, macroturbulence, pressure broadening, and 
eclipse blockage on the model absorption line kernel. 

The kernel is then convolved with a line-list in the 
wavelength region between 4460 A and 4485 A. The po- 
sitions and line strengths are obtained from the Vienna 
Atomic Line Database (VALDlj for t he stellar param- 
eters given by Holmgren et al.l TlQQO) and are listed in 

" Tmir 



Table [2| The He I, Mgll, Allirana some of the O II 
lines consist of multiple lines with energy levels spaced 
closely enough (< 0.3 A) that we model them as single 
lines. We also omitted one S III line which has nearly the 
same wavelength as one of the O II lines. 



^ http : //ams . astro . univie . ac . at/vald/ 
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Fig. 6. — Observations of NY Cep outside of eclipses. These panels show spectra of both stars in the NYCep system in the 
spectral region around the He I line at 4471 A and Mgll line at 4481 A. The number in each panel indicates NYCep's orbital phase at time 
of observation, defined such that phase corresponds to periastron. The gray solid lines represent the data and the (red) dashed and (blue) 
dash-dotted lines are the simulated absorption lines of the primary and secondary, respectively. The black line is the best fitting model. 
The dots around the line at a fiux level of 1.05 represent the differences between the data and the model. 



2.2. Eclipse timing 
Since we did not obtain new photometri c data, we re 



Ahnl ([1992]) (in their Ta- 
period and time of pri- 



hed on the data presented by 
ble 2) to constrain the orbita" 
mary minimum hght. The lack of error estimates by |Ahn| 
|l992 ) presented a comphcation, which we dealt with as 
follows. We fitted a linear function of epoch to the re- 
ported times, assuming equal errors in all the measure- 
ments, determined by the requirement = ^dof- The 
resulting error in each time was 0.012 days, and the re- 
sulting ephemeris are given in Table [ 3 } Our re sult for the 
period agrees with that given by lAhnl (|1992D, although 
we found a slightly different time of primary minimum 
light, presumably due to different weighting of the mea- 
surements. 



TABLE 3 

Results eclipse timings 



Parameter 



This work 



"Ahn|(|l992l) 



Period 15^27566±0.00002 15^27566 
Tmin,i 1973 [HJD-2 400 000] 41903.819±0.006 41903.8161 
Tminj 2009 [H JD-2 400 OOP] 55117.265^0.016 



Then we used a Monte Carlo approach to calculate the 
predicted times of minimum light at the epochs of our 
observations, in a manner that respects the correlation 
between the parameters of the ephemeris. We created 
10^ fake data sets by ad ding Gaus sian perturbations to 
the times of minima of | Ahn (1992), with a standard de- 
viation of 0.012 days. We then fitted a line to each of 



these fake data sets and used the linear fit to calculate 
the expected mideclipse time on 2009 Oct. 12/13, the 
last night of our eclipse observations. The mean and 
standard deviation among all 10^ results were taken to 
be the value and the error in the mideclipse time. The 
results are given in Table [3| 

2.3. The fitting procedure 

We now describe the model parameters in detail. Seven 
parameters describe the Keplerian orbit: the time of min- 
imum light (Tmin,!), the eccentricity (e), the argument of 
periastron (cj), the velocity semiamplitudes of the pri- 
mary and secondary (Kp and Kg), and the velocity off- 
sets (7p and 7s) jj We held the period (P) fixed, as its 
uncertainty is negligible over the 3 months spanned by 
our observations. 

Another 4 free parameters specify the photometric as- 
pects of the eclipse: the light ratio between the two stars 
at the wavelength of interest (Lg/Lp at 4500 A), the frac- 
tional radii of the stars (rp and rg), and the orbital in- 
clination {io). The linear limb darkening parameter (ui) 
was held fixed at 0.4 for both stars (Gray 2005). 

The model of the velocity fields introduces 6 param- 
eters: the projected equatorial rotation speeds ('usinip 

^ The velocity offsets are approximately equal to the radial ve- 
locity of the NYCep barycenter with respect to the Solar system 
barycenter, but the offsets can differ from each other due to subtle 
factors specific to each star, such as the gravitational redshift and 
line blending. The latter term refers to the fact that the stars have 
slightly different strengths and positions of closely spaced lines (un- 
resolved in our spectra), due to the differences in effective temper- 
ature, surface gravity, and other atmospheric parameters. These 
differences can result in a systematic bias in the inferred systemic 
velocity. 
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Fig. 7. — Observations of NY Cep during primary eclipse. Similar to Figure [g] but for spectra obtained during primary eclipse. 
The number in each panel indicates the time from mideclipse, in hours. The corner of each panel shows an illustration of the eclipse phase. 
In the best-fitting model, /3p = 2 ±4°, indicating a close spin-orbit alignment. Focusing on the Mgll line, one can see in particular for the 
observations at — 2.8h and 2.8 h how the RM-effect creates a redshift during the first half of the eclipse and a blue shift during the second 
half of the eclipse. For the observation at — 0.3h, near mideclipse, blue and redshifted light are blocked by similar amounts. 



and V sin is), the Gaussian width of the macroturbu- 
lence for the primary star {CrtP)^ the half- width at half- 
maximum of the Lorentizan function representing pres- 
sure broadening for each star (^p and ^g), and finally, the 
parameters of greatest interest for this study, the sky- 
projected spin-orbit angle ( 0r^). The a n^le is defined ac- 
cording to the convention of Hosokawa (1953), such that 
/3 = 0° when the axes are parallel, = ±90" when they 
are perpendicular, and /3 > when the time-integrated 
RM effect is a redshift for an orbit with io < 90° The 
macroturbulence parameter for the secondary was held 
fixed to 20kms~^, since the secondary is significantly 
fainter than the primary, causing the results to be insen- 
sitive to this parameterjj 

Eight additional parameters are needed for each star 
to describe the relative depth of the spectral lines in the 
wavelength range being modeled. In addition, 3 free pa- 
rameters are needed to describe the quadratic function 
used to normalize the continuum level of each of the 46 
spectra. The total number of parameters is therefore 
7 (orbital) + 4 (photometric) + 6 (velocity fields) + 
8x2 (line depths) + 3 x 46 (normalization), for a to- 
tal of 171 adjustable parameters. Although this may 
seem like a large number, the 138 normalization param- 
eters are "trivial" in the sense that they can be opti- 



^ The parameter A, which was introduced bylOhta et al.|(|2005fc 
and is commonly used in the exoplanet community, is simply —f:^. 

^ Indeed, although we included macroturbulence for complete- 
ness, the key results for /3p do not depend on the details, and 
similar results were obtained even when macroturbulence was ig- 
nored. One can see from Figure^ that macroturbulence has only 
a small effect on the overall shape of the model absorption line, 
given the high v sin i of the primary. This is even more true of the 
secondary, which rotates about twice as fast as the primary. 



mized separately in subsets of 3, for a given choice of the 
other parameters. Also, several of the other parameters 
are subject to a priori constraints. More detail on these 
points is given below. 
The fitting statistic was 



46 213 



i=i j=i 



fx,^j{ohs) - /A,,.(calc) 



T^ini - 2455117.265 y 



78 



0.016 

2 



Ls/Lp - 0.37 
0.06 



0.7 



0.5 



where /a • ■ (obs) is the observed fiux in pixel j in obser- 
vation i and fx. . (calc) is the calculated fiux based on a 
particular choice of the model parameters, and (Ji rep- 
resents the uncertainty in the fiux pixel for observation 
i. The latter was taken to be the reciprocal of the SNR, 
estimated as described in section [2] In this equation, the 
first two terms represent the usual sum-of-squares, and 
the remaining t erms represent Gauss ian priors based on 
the results of Holmgren et al.' ("1990) (except for Tmin,i, 
which was discussed in the previous section). 

To optimize the parameters we used Levenberg- 
Marquardt (LM) least-squares minimization. The main 
optimization was conducted for the nontrivial parame- 
ters. However, whenever was computed within that 
process, the values of the trivial (normalization) param- 
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Fig. 9. — Apparent radial velocities of NY Cep. a: The 

apparent radial velocity (RV) of the primary (red filled circles) and 
secondary (blue open circles) as a function of orbital phase. The 
line positions were determined by fitting Gaussian functions. The 
solid line is the calculated radial velocity based on our model to 
the shape ob the absorption lines, including Keplerian motion and 
the RM effect. b,c: Close-ups of the RM effect during the primary 
eclipse. The upper panel shows the RM effect plus orbital motion 
and the lower panel shows only the RM effect, a redshift during the 
first half of the eclipse, and a blueshift during the second half. One 
can see that a fit to the RV data would overestimate the amplitude 
of RM signal. We emphasize that the eclipse radial velocities are 
shown here for illustration only, as a concise visual summary of the 
complex distortions of the eclipse spectra: they were not used in 
our quantitative calculations. 



Fig. 8. — Observations of NY Cep during primary eclipse. 

Top. — Grayscale depiction of the spectra obtained during and after 
primary eclipse. Each spectrum is shifted into the rest frame of the 
primary star, and the best-fitting model of the secondary spectrum 
has been subtracted. The most prominent lines are He I (seen here 
blended with the weaker O II lines) and Mg II 4481 A . Dashed lines 
indicate the beginning and end of the eclipse. Middle. — The 
same as the top panel, after subtracting the mean absorption line 
profile shown in Figure [3] The RM effect is evident as the white 
residual traveling from the blue to the red as the eclipse progresses. 
Bottom. — Residuals between best fitting model and data. 



eters were first optimized using a separate 3-parameter 
minimization for each order. This process is equivalent 
to the "Hyperplane Least Squar es" method that was 
described and tested by ,Bakos et al.| (2010). To esti- 
mate the parameter uncertainties, we us ed the bootstrap 



Press et al. 



(1992), with 5x10^ re- 



method described by 
alizations. 

With 171 parameters, 5 Gaussian priors, and 46 x 213 
data pixels, there are 9632 effective degrees of freedom. 
The best-fitting model has = 9148, or = 0.95. 
The low is probably the result of underestimating the 
SNR of the relevant lines. The estimate was based on 
the continuum nearer to the edges of the order, where 
the SNR is lower (see Figure |2|. We have not attempted 
to correct for this effect. 

3. RESULTS 

The results for the model parameters are given in Ta- 
ble [4j Figure [6] shows the model fitted to a few out- 
of-eclipse spectra in the vicinity of the He I and Mgll 
lines. Figure [7| shows the same for a few of the spectra 
obtained during primary eclipse. A grayscale represen- 
tation of the eclipse spectra and post-egress spectra is 
shown in figure [8| We also display the apparent radial 



velocities in NY Cep in figure [9| 

3.1. Orbital parameters 

Our result for Tmini is 1.2cr away from the value calcu- 
lated from the older eclipse timings. Since the ephemeris 
is based on observations obtained more than 20 years 
ago, and we do not know the true uncertainty in each 
measurement, we regard this as good agreement. Like- 
wise the eccent ricity is within 2a of the value given by 
'Holm gren et al.| ( |1990[ ) . The measured value for uj agrees 
with the previous value, but as it is expected to change 
over time due to apsidal motion a more careful compar- 
ison is needed, as discussed in section [4j 

Our results for Kp and 7 n agree with the values found 
by Holmgren et al.' (1990). However, our results for the 
secondary star are different: we find a lower value for Kg, 
and we fou nd A7 = 7^ — 7 n to be consistent with zero. In 
contrast, iHolmgren et al. (1990) found the secondary to 



be redshifted by 26km s~^ relative to the primary, and 
Heard & Fernie (1968) found an even greater difference of 
32 km s~^ . We speculate that the previous measurements 
were subject to bias because of the relative faintness of 
the secondary, its rapid rotation, and the need to use 
pressure-broadened (and asymmetric) lines. Based on 
the agreement in 7p it seems that the systemic velocity 
has not changed significantly over the last decades. 

3.2. Stellar parameters 

We further find with our fit to all spectra a luminosity 
ratio (Ls/Lp) of 0.36 ± 0.05 at 4500A. This result agrees 
with the value 0.35 that is obtained from comparing the 
equivalent widths of the He I (4471 A) line profiles shown 
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TABLE 4 

Parameters of the NY Cephei system. 



Parameter 


This work 


Literature values 


Orbital parameters 


Time of primary minimum, Tmin,i [HJD-2 400 000] 


55117.287±0.009 


41903.816lt 


Period, P [days] 


15.27566 (fixed) 


15.27566t 


Eccentricity, e 


0.443±0.005 


0.48±0.02* 


Argument of periastron, cu [deg] 


56.3±1 


58±2* 


Orbital inclination, io [deg] 


78.8±0.7 


78±1* 


Velocity semiamplitude (primary), Kp [kms"-*^] 


113.8±1.2 


112±2* 


Velocity semiamplitude (secondary). Kg [kms"-*^] 


139±4 


158±8* 


Velocity offset (primary), 7p [kms"-*^] 


-17.3±0.7 


-15±2* 


Velocity offset (secondary), 7s [kms"-*^] 


-21±4 


9±6* 



Stellar parameters 



Light ratio, Lg/Lp 


0.36±0.05 


0.37±0.06* 


Fractional radius (primary), rp 


0.086±0.007 




Fractional radius (secondary), rg 


0.084±0.009 




Projected rotation speed (primary), -usinip [kms"-*^] 


78±3 


75±10* 


Projected rotation speed (secondary), i;sinzs [kms"-*^] 


155±6 


125±14* 


Macroturbulence parameter (primary), (^p [kms"-*^] 


23±6 




Macroturbulence parameter (secondary), [kms"-*^] 


20 (fixed) 




Pressure broadening parameter (primary), ^p [kms"-*^] 


36±1 




Pressure broadening parameter (secondary), [kms"-*^] 


54±6 




Linear limb darkening parameter (primary), Up 


0.4 (fixed) 




Linear limb darkening parameter (secondary), Us 


0.4 (fixed) 




Projected spin-orbit angle (primary), /3p [°] 


2±4 





Indirectly derived parameters 



Projected orbital semimajor axis, a sin 

Mp sin^ i [Mq] 

Ms sin^ i [Mq] 

Primary mass, Mp [Mq] 

Secondary mass, Mg [Mq] 

Primary radius, Rp [Rq] 

Secondary radius, Rs [Rq] 



[Rq] 



68.6±1 
10 1+°-^ 



8.3^ 



^+0.5 
-0.3 

10.7±0-^ 

8.81^4 
6.0±0.5 
5.8±0.5 



71.4±2.2* 
12±1* 
8.7±0.6* 
13±1* 
9±1* 
6.8±0.7* 
5.4±0.5* 



t Data from 
i<: Data from 



Ahn|(|l992|> 



Holmgren et al. ( 1990 1 



Both values are consistent with 0.37 ± 0.06 



olmgren et al. (1990), based on spectropho- 



in Figure 3^ 
found bylHc 
tometry. 

Our model gives projected rotation speeds of 78 ± 1 
kms~^ for the primary 155 ± 4 kms~^ for the sec- 
ondary. The width of the Gaussian function describing 
the macroturbulence in the primary is 23 ± 3 kms~^, 
and the half-width at half-maximum of the Lorentzians 
describing pressure broadening are 36 ± 1 kms~^ and 
54 ± 6 kms-^ These error estimates must be under- 
stood as internal to our model, which we recognize may 
not be completely realistic, especially as it pertains to 
pressure broadening. We have assumed that absorption 
lines form at a specific value of pressure, which is not nec- 
essarily true. For this reason, the Lorentzian widths have 
no simple physical interpretation, and the results may be 
somewhat biased for all other parameters affecting line 
broadening (namely vs'mi and (). For comparisons to 



^ In estimating the luminosity ratio from the EWs, we accounted 
for the small increase (~10%) in the EW of this line towards lower 
temperatures in the spectral range from to BO to B2. See Figure 1 
of Leone & Lanzafame (1998). 



other analyses obtained with different instruments and 
different analysis procedures, we recommend using the 
more conservative error estimates given in Table [4j For 
the projected rotati on rates, we find a sim ilar result for 
the primary as did Holmgren et al. (1990), but we find a 
significantly higher v sin i for the secondary, perhaps for 



the same reasons given in section 3.1 



As for the projected spin-orbit angle /3p, the focus of 
this study, we find /3p = 2 ± 4° . There is no correlation 
between this parameter and the linewidth parameters, 
so the concerns raised above about the oversimplified 
pressure-broadening model do not apply here. A strong 
correlation does exist with the time of minimum light 
(see Figure 10). This implies that future photometric 
observations to refine the eclipse ephemeris would lead 
to smaller uncertainty in /3p, although the uncertainty is 
already quite small. 

The absolute radii of the stars can be calculated from 
the fractional radii and the other orbital parameters. We 
find radii consistent with the values given in the litera- 
ture, Rp = 6.O±O.5i?0 and = 5.8±O.6i?0, which is 
not surprising since our results were strongly infiuenced 
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Fig. 10. — Results for time of minimum light and /3p, based 
on our analysis of the stellar spectra using the bootstrap method. 
There is a strong correlation between these two parameters. 



by the priors on those parameters. For the stellar masses, 
we find Mp = 10.7±0.8 M© and Ms = 8.8±0.5 M© for the 
primary and secondary, respectively. The value found for 
the primary is smaller than that reported previously, due 
to the previously mentioned discrepancy in the velocity 
semiamplitude of the secondary. 

4. DISCUSSION 

Primary spin axis — We have measured the angle be- 
tween the projections of the orbital and stellar angular 
momentum vectors for the primary star in NYCep and 
find that these projections are aligned within the uncer- 
tainty of our measurement (/3p = 2 zb 4°). We believe 
this result to be robust, even if our model is simplified 
in many respects. We have neglected any changes in the 
limb- darkening law within the lines as compared to the 
continuum, as well as a number of other "second-order" 
effects, such as gravity brightening and differential rota- 
tion. We experimented with more complex models in- 
cluding these phenomena and found that they produced 
effects too small to affect our conclusions about spin- 
orbit alignment. Likewise we experimented with a "two- 
layer" mo del for the pressure-broadened lines [see fig- 
ure 11.5 of |Gra y|(r2QQ5)] and found that while it provided 
a slightly better fit to the data, it introduced several new 
free parameters and led to no significant changes in any 
other parameters. 

Our result for /3p gives a lower bound on the true angle 
between the stellar and orbital spins, the obliquity (i/j). 
We however expect the obliquity to be not much larger 
than P unless the primary spin axis is highly inclined 
towards the line of sight, an unlikely scenario. 

Absolute dimensions — For our di scussion we ado pt the 
effective temperatures derived by Holmgren et al. ( 1990) , 
which are printed in table [2] Because we do not have 
any new photometry, we do not attempt to derive new 
effective temperatures or a new age for the system. We 
adopt the previous classification of NY Cep as a young 
system (~6 Myr) at a distance of 750 ± 90 pc, making 
it consistent with membership in the Cepheus OB III 
association. Our finding of a lower mass for the primary 
does not change the picture of NY Cep significantly given 
the uncertainties in other parameters. 



Stellar rotation — If the spin periods were equal to the 
orbital period, then the stellar rotation velocities would 
be 20 kms~^ and 19 kms~^, for the primary and sec- 
ondary, respectively. These numbers are smaller than 
our results for vsinz, which are 78±3kms~^ and 154 ±6 
kms~^. Assuming = for both stars, the implied 
rotation periods are 3.9 and 1.9 days. 

The stars do not seem to be pseudosynchronized, ei- 
ther. The ratio of rotational to orbital frequencies is 3.9 
and 8.0 for the primary and secondary, respectively. In 
the |Hut| (p'981) model of pseudosynchronization, the ra- 
tio would be 2.3, which is lower than the observed ratios. 
If the stars had been significantly infiuenced by tidal evo- 
lution, then one would expect the rotational frequency to 
be smaller for the secondary star than for the primary 
star, because the less massive star should synchronize 
first. In contradiction with this expectation, we observe 
the secondary to have a higher rotation frequency than 
the primary. 

All of these findings are in agreement with a picture in 
which the system is only a few million years old, and tidal 
forces have not had enough time to alter the rotational 
state of the stars after their arrival on the zero-age main 
sequence. 

Apsidal motion — With the data at hand it is not possible 
to investigate the internal structures of the two stars in a 
stringent way via observations of apsidal motion. Firstly, 
the uncertainties in the radii and ages are so high that 
the theoretical expected apsidal motion has a high un- 
certainty. Using the system parameters, together with 
internal structure constants of k2 = —0.080 ± 0.3^ for 
both stars, we derive an expected apsidal motion rate 
of 19^9^ arcsec per cycle. Secondly, the measurements 
of the argument of the periastron (uj) are also quite un- 
certain. Using the values from |Heard fc Fernie| (|1968[), 
who derived an u of 50.9 ± 3.7, from [Holmgren et al. 



(1990), and from this work, we derive a measured apsi- 
dal motion rate of 6 ± 11 arcsec per cycle. Therefore the 
expected and measured apsidal motions are consistent, 
but the bounds are too weak to give additional informa- 
tion. To make progress on this field both the expected 
and measured values for the apsidal motion have to have 
greater precision and accuracy. 

5. CONCLUSION 

Through an analysis of high resolution spectra of the 
NY Cep system, observed during and outside of primary 
eclipse, we have calculated the orbital and primary and 
secondary parameters. We also reanalyzed times of min- 
ima from the literature to obtain an updated ephemeris 
for the epoch of our spectroscopic observations. We find 
a ~15% smaller mass for the primary star than earlier 
studies, and derive results for the other parameters that 
are similar to those reported in the literature values. 

We measured the angle on the sky between the primary 
spin axis and the orbital axis by exploiting the Rossiter- 
McLaughlin (RM) effect. We find that the projections of 
these axes on the sky are closely aligned (/3p = 2±4). The 
close alignment does not seem to have been the result of 
tidal evolution, because tidal evolution would also have 

^ We estimated the apsidal motion constants using tables from 
|Claret|fl995 ). 
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synchronized the rotational and orbital periods, which is 
not observed. 

The finding of a close spin-orbit alignment is in strong 
contrast to the situation in the other young, early type, 
detached close binary, DIHer, for which the orientation 
of the axes has been measured (Paper II). For that sys- 
tem we found strongly misaligned spin and orbital axes. 
While NYCep and DIHer are similar in that they both 
harbor young, well detached, early type stars, and have 
high orbital eccentricities their formation and/or evolu- 
tion seem to have taken quiet different routes. 

DIHer and NYCep are currently the only early-type 
detached binaries for which the orientation of the stellar 
spin axes is known. Both systems are on eccentric or- 
bits, have detached components, and are young enough 
to exclude any development of the spin-orbit alignment 
during main-sequence life due to tidal interactions. Yet 
the axes in the DIHer system are strongly misaligned 
(/3p = 72 ± 4; /3s = —84 ± 8) while the primary axis in 
NY Cep is well aligned (/3p = 2±4). Indeed, the projected 
alignment of NY Cep is even closer than the alignment 
between the Sun's spin axis and the ecliptic plane (7°). 

As of yet there is no good explanation for the misalign- 
ment of DIHer. Broadly speaking, it could be primor- 
dial, or it could be the result of an interaction between 
the stars or between the stars and the disk from which 
they formed, or it could be due to torques from addi- 
tional bodies in the system. The results in this paper do 
not point to a particular mechanism, but they do mean 
that any successful theory must pass the "NYCep" test, 
namely, it must explain the difference between those two 



systems. 

Obtaining measurements of spin orbit angles in a num- 
ber of early type close binary systems might reveal which 
environmental variables determine the orientation of stel- 
lar rotation axes in these systems. For example if the 
interaction with a distant companion via the Kozai mi- 



gration (Fabrycky & Tremaine||2007 ) formed close bina- 



ries, than misaligned systems should be more likely to 
have a detected companion and closer systems might be 
more likely to have misaligned axes than wider systems. 
If alignment is a simple function of coherence length dur- 
ing star formation, then the opposite should be true: 
the closer the system, the stronger the tendency toward 
alignment. 

This is why we have started a small survey to measure 
the orientation of stellar spin axes in other early-type 
binaries. 
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